Counterion desorption and globule-coil transition of a polyelectrolyte under tension 
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We explore the mechanical response of a single polyelectrolyte chain under tension in good and 
poor solvents using a combination of simulation and theory. In poor solvents, where the equilibrium 
state of the chain is a collapsed globule, we find that the chain undergoes a globule-coil transition, as 
the magnitude of the force is increased beyond a critical value. This transition, where the polymer 
size changes discontinuously from a small to a large value, is accompanied by release of bound 
counterions from the chain. We explain these results by adhering to a statistical mechanical theory 
of counter-ion condensation on flexible polyelectrolytes. 
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A polyelectrolyte (PE) is a charged polymer with ion- 
izable side groups on its backbone that dissociate in polar 
solvents, resulting in a macro-ion and oppositely charged 
counterions In solution, the size of a polyelectrolyte 
is intricately coupled to its effective charge. The conse- 
quence of this coupling is most significant when a poly- 
electrolyte undergoes a globule-coil transition. Such a 
phase transition can be effected by either changing the 
solvent quality (temperature) or applying a tensile force. 
In this letter, we employ Langevin dynamics simulations 
to examine the elastic response of a flexible polyelec- 
trolyte chain that undergoes a globule-coil transition un- 
der externally applied tensile forces. We explain the sim- 
ulation results using a theoretical calculation based on a 
double minimization of the PE free energy [2|. 

Our work is motivated by single molecule force spec- 
troscopy (SMFS) experiments on biologically occurring 
polyelectrolytes [3[, such as proteins and DNA, where 
AFM and optical tweezers are used to examine their 
elastic response. The elastic response obtained from 
these experiments provides insight into their consti- 
tutive behavior bonding interactions and bind- 
ing energies [1]. Although these experiments can exam- 
ine conformational changes of biopolymers, the role of 
charges in such transitions cannot be experimentally de- 
termined. However, using simulations^ — 0, we can 
track polymer conformations and calculate its degree of 
ionization simultaneously. 

In a poor solvent, a polyelectrolyte (PE) collapses 
into a globular configuration due to unfavorable solvent- 
polymer interactions. However, in order to minimize elec- 
trostatic repulsion between monomers, most of the coun- 
terions are adsorbed on the chain contour (Fig. [ija)). In 
a good solvent, the PE adopts a swollen coil-like con- 
figuration because of favorable solvent-polymer interac- 
tions, and repulsive electrostatic interactions between 
monomers. The counterions now have the freedom to 
maximize their translational entropy, and some desorb 
from the chain to become mobile in solution. The fraction 
of mobile counterions, with respect to those adsorbed on 
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FIG. 1: Configuration snapshots of a single polyelectrolyte 
chain (red) under tension and counterions (blue). Panel (a) 
shows a PE globule with adsorbed counterions. The globular 
configuration just prior to the transition is shown in panel 

(b) while the coil configuration beyond a critical value of the 
applied force is shown in panel (c). Panel (d) shows a PE 
chain stretched to near complete extension at a higher value 
of applied force. 

the chain, depends on the relative strengths of electro- 
static interactions and thermal fiuctuations. 

We simulate the stretching of a PE chain under equal 
and opposite tensile forces applied to its ends. Our main 
result is that at a critical value of applied force the PE 
chain undergoes a globule-coil transition (Fig[T](b) and 

(c) ), characterized by a discontinuous change in chain 
size. This transition is accompanied by a correspond- 
ing discontinuous increase in the degree of ionization (to- 
tal number of free counterions in solution) of the chain. 
This is the first demonstration of counterion desorption 
accompanying a globule-coil transition, and calls for a 
careful analysis of previous theories where the ionization 
of the chain is assumed to be a parameter that remains 
constant throughout the transition [loj . We explain our 
simulation results using a variational formulation of the 
PE free energy that incorporates the stretching of the 
polymer under tension. Comparing the relative strengths 
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of the different terms in the PE free energy, we find that 
at the transition, gain in mobile counterion entropy is 
greater than the chain-counterion adsorption energy and 
leads to counterion desorption. 

The polyelectrolyte is modeled as a freely-jointed bead- 
spring chain comprised of = 100 spherical beads con- 
nected by iV — 1 springs, where the length of each spring 
corresponds to a Kuhn step, Iq. The total charge Q of 
the polyelectrolyte is uniformly distributed among all the 
beads (charge on each bead, q = Q/N). Charge neutral- 
ity is maintained by introducing N monovalent counte- 
rions, modeled as spheres, each with a point unit elec- 
tric charge —q. The polyelectrolyte and counterions are 
placed in a medium of uniform dielectric constant. The 
system is enclosed in a cubic box of edge length 50/o 
with periodic boundary conditions applied in each coor- 
dinate direction. Each bead represents a Brownian parti- 
cle, whose motion is governed by the Langevin equation, 
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where rrii, and represent the mass, position and ve- 
locity of particle z, respectively. The random force F[(t) 
acting on each bead satisfies the fluctuation-dissipation 
theorem, (F[(t)FJ(t')) = &CkBT5{t ~ t')5,j, where ks 
is the Boltzmann constant, T the absolute temperature, 
and C is the frictional drag. Externally applied forces, 
Ff^*, of equal magnitude, but opposite directions act 
on the terminal beads of the chain. However, Ff^* = 
for the remaining polymer beads. Equation [1] was solved 
numerically using the velocity- Verlet algorithm, with a 
time step At = O.OOIt, where t = lo[m/ eLjY^'^ ■ 
The net interaction potential, U{rij), is given by. 



U{nj) = UlJ + Ubond + Uelec, 



(2) 



corresponding to excluded volume, bond stretching and 
electrostatic interactions, respectively. The excluded vol- 
ume interaction between any two beads is described by. 
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where Vij = jr; — rj|, e^j is the Lennard- Jones interaction 
parameter and a is the bead diameter. The bead diam- 
eters, (T, used in Eq. ([3]) for bead-bead, bead-counterion, 
and counterion-counterion interactions are 0.75 Iq, 0.55 
Iq and 0.35 Iq, respectively [lH. We use e^j, lo and m 
as scales for energy, length, and mass, respectively. The 
masses of individual counterions and polymer beads are 
set to one half and a unit mass, respectively The bond- 
ing interaction between adjacent beads is governed by 
the FENE potential given by. 



Ubondin) 



1 



-krh In 



1 - {n/RmaxY 



(4) 






Globule 


Coil 





0.75 

e 



u 



FIG. 2: The variation of radius of gyration of the polyelec- 
trolyte chain Rq as a function of the solvent quality elj show- 
ing a coil to globule transition at eL.j ~ 0.6. The F values were 
varied between 0.1 and 7.0, corresponding to low and high 
electrostatic strengths, respectively. Similarly, e_Lj was varied 
between 0.1 and 2.5 corresponding to good and poor solvents, 
respectively. We compute the size of the chain, R^/N, and 
its degree of ionization, a, for these parameter ranges. 



where is the separation distance between adjacent 
beads, Rmax — 1.5Zo is the maximum allowable separa- 
tion distance between bonded beads and k is the spring 
constant. The electrostatic interactions between beads is 
given by the Coulombic potential 



Uelec ij^ij ) 



ATTeoerTi 



(5) 



where eq is the vacuum permittivity and is the rel- 
ative dielectric constant of the medium. The particle- 
particle particle-mesh (PPPM) method with a precision 
of 10^'^ was utilized to calculate the long-range electro- 
static interactions in the system. Electrostatic interac- 
tions are parameterized by the Coulombic strength pa- 
rameter F = t^-, where In = i — — ^ is the Bierrum 
length. The simulations were executed in the open-source 
molecular dynamics simulation package LAMMPSfl^. 

Shown in Fig. [2] is the variation of the normalized ra- 
dius of gyration, (Rg) /N, of the chain as a function of the 
solvent quality, e^j, at F = 4.0 (for bio-polymers, such 
as DNA, F hes between 2.0 and 4.o[ii|). A coil to glob- 
ule transition, signalled by a large change in {R'^)/N is 
observed with increasing e^j (or decreasing solvent qual- 
ity), with the transition occurring around e^j w 0.6. We 
note that for poor solvents clj ^ 0.75, the chain collapses 
into a globule of size smaller than Gaussian dimension. 
Inset shows the phase diagram in the F — e^j plane, de- 
marcating the globule and coil phases. The phase bound- 
ary is obtained by calculating (Rg) /N as a function of 
these parameters. 

In order to extract the mechanical response of a sin- 
gle PE chain, it is subjected to a pulling and retraction 
force cycle, and its root-mean-square end-to-end distance 
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FIG. 3: Force vs. extension (top panel) and vs. degree of ion- 
ization (bottom panel), for F = 4.0, and as the solvent quality 
is varied during pulling (circles) and retraction (triangles). At 
zero force the chain is in a globular state for eLj ~ 0.9, and 
1.5, while it is in the coil phase for el./ = 0.1, and 0.5. The 
degree of ionization jumps from a ~ 0.12 to a ~ 0.23 at the 
coil-globule transition transition at fc ~ 3.0 (pulling cycle) 
indicating counterion desorption. 



FIG. 4: Force vs. extension (a) and vs. degree of ioniza- 
tion (b) from variational calculation. Depending on tem- 
perature a globule- rod transition (T = 168, solid), globule- 
coil transition (T = 204, dotted) and a coil-rod transition 
(T = 254, dashed) is seen. Top inset shows the total free 
energy vs. extension, with minima corresponding to glob- 
ule (/ = 0,ri = 0.0002, a = 0.0151) (top) and coil (/ = 
0.0003, h = 33, a = 0.19) (bottom) states at T = 252 K. 



is monitored as a function of time. The force is incremen- 
tally increased from / = to / = 10 (/ = F/Fq, where 
Fo = kBT/lo, sets the scale of force) in the pulling cycle, 
and subsequently brought back to / = in the retrac- 
tion cycle. For a given value of applied force, the PE 
chain is equilibrated for 5 x 10^ time steps. The normal- 
ized chain extension, (x) /NIq , is the time averaged root- 
mean-square end-to-end distance calculated over this du- 
ration. Simultaneously, the degree of ionization of the PE 
chain, defined as, 
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where (ric) is the average number of counterions con- 
densed on the polymer ll|, is also calculated over the 
force cycle. 

Figure [3] shows the extension of the chain, {x)/NIq, 
(Fig.[3{a)), and its degree of ionization, a, (Fig.[3]Jb)), as 
a function of externally applied force, / at T = 4.0 and 
ELJ = 0.1,0.5,0.9,1.5. For clj = 0.1,0.5, corresponding 
to a good solvent, (x) /NIq increases continuously from a 
swollen coil configuration to its maximum extent, as / in- 
creases from 0.1 to 10. In both cases, the force-extension 
curves for pulling and retraction are identical. In con- 
trast, in a poor solvent (elj = 1.5), the chain exists as 
a compact globule at / = (sec Fig. Hlja)). As / is in- 
creased, (x) /NIq remains nearly constant until the force 
reaches a threshold value, fc ~ 3.0. Beyond this value 
a large increase in {x) /NIq is observed for a very small 



increment in applied force (Figs. [Ijb) and (c)), defining 
the force-induced globule-to-coil transition for a PE in a 
poor solvent. 

A similar trend is observed in the variation of the de- 
gree of ionization of the PE (see Fig. E^b)) with the ap- 
plied tensile force. For e^j ~ 1.5, a also shows a jump 
discontinuity at the same value of critical force, fc ~ 3.0. 
The degree of ionization, a, approximately doubles from 
a value of 0.12 to 0.23 at the transition. This suggests 
that the dramatic increase in the size of the chain at the 
transition leads to a reduced electrostatic repulsion be- 
tween polymer segments. It is therefore favorable for the 
counterions to unbind from the chain and become mo- 
bile, thus gaining translational entropy that leads to a 
lowering of the system free energy. 

We also note that both the force-extension curves show 
hysteresis for e^j = 1.5. This is a common feature of first 
order phase transitions where the system gets trapped in 
a metastable state, whose lifetime is longer than the ob- 
servation time for simulations. The degree of ionization 
as a function of force shows similar hysteresis, indicat- 
ing that the counterions relax on time scales much faster 
compared to the polymer relaxation time. Distinguish- 
ing long-lived metastable states from equilibrium phases 
is difficult, and requires a calculation of the free energies 
of the different phases. 

A variational formulation of the free energy of a PE has 
been recently proposedH and analyzed in the context of 
globule-coil transition[13[. The total free energy of the 
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PE system is given by, 

Ftot = Fcion + FpE + Fads + Ften ■ 



(7) 



Counterion contribution to the free energy, represented 
by the first term, 

Fc/NksT = a log a + (1 - a) log(l - a) 

+ {{ap + Cs) log(ap + Cs) + Cs logc^ - {ap + 2cs))/p 



ap 



2c, 



n3/2 



(8) 



where a is the degree of ionization of the chain, the 
salt concentration, p the monomer density, and Ib the 
Bjerrum length normalized by the Kuhn segment length 
Ib — Ib/Io- The first two terms in Eql8]correspond to the 
translational entropy of the adsorbed and mobile counte- 
rions respectively, while the third term is due to fluctua- 
tion of free ion density as given by Debye-Hiickel theory. 
The free energy of the polyelectrolyte, obtained from the 
Edwards Hamiltonian[^ by a variational calculation is 
given by 
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where (i?^) = Nl^h, and a = k'^NIi/6, k being 
the inverse Debye screening length given by = 
AiTlB{ap + 2cs). The last term in Eq. [5] arises 
from the electrostatic interactions between segments 
of the polyelectrolyte chain and is given by 6o(a) = 
^ - ^) exp[a]erfc[yS] + i+ - -V| _ 
while the first three terms correspond to configurational 
entropy of the chain, a repulsive excluded volume and 
steric three body interaction terms respectively. 

Electrostatic interaction energy of ion-pairs formed due 
to counterions adsorbed on the chain is given by 



Fads/NkBT^-{l-a)8h 



(10) 



where 5 = {e/ei){lo/d,), is the ratio of the bulk dielectric 
constant e, to ei the local dielectric constant and d is 
the length of the dipole. Contribution to the total free 
energy due to interactions between ion-pairs marginally 
affect the coil globule transition and are hence ignored. 

In addition to these contributions the tensile force ap- 
plied to the ends of the polyelectrolyte chain has a con- 
tribution. 



Ften/NkBT : 




(11) 



to the PE free energy 1J| . We obtain the chain extension 
ratio li, and the degree of ionization a as functions of 



force / by a double minimisation of the total free energy 
Ftot- 

Fig. 3] summarizes our theoretical results. We consider 
a PE of iV = 1000, at Cs = 0, p = 0.0005 and S = 3.0. 
The excluded volume interaction parameter is taken as, 
w = 1 — 6/T. As an example, we consider 6 = 400/^ and 
Ib = 900/T. For the parameters used, a discontinuous 
globule-coil transition is observed for W3 < = 8.5[l3j. 
In Fig. IHa), we plot li as a function of / for = 0.05 
and T = 168, 204, 254 K. At T = 254 K, the equihbrium 
coil state continuously stretches to a rod (^i = TV, used 
as an upper bound in our theory) with increasing /. At 
a lower temperature, T ~ 204 K, the globular PE under- 
goes a discontinuous globule-coil transition at a critical 
force, /i^' « 0.002. On further increase in /, the coil 
stretches to a rod. A different discontinuous transition is 

observed for T = 168 K, where the globule transitions to 
(2) 

a rod at a critical fc ~ 0.003 bypassing the coil phase. 
This is because the coil phase is thermodynamically un- 
stable for these temperatures [l3| . 

The PE free energy (at T = 252 K) as a function of h is 
shown in Fig.jH^a) inset. At / = 0, the global free energy 
minimum occurs at h « 0.002, corresponding to a globu- 
lar PE (inset top panel). For f > fc = 0.0003 the global 
minimum in free energy shifts to li w 33 corresponding to 
a coil phase. At this force, the limit of globule metasta- 
bility is reached, leading to a discontinuous first-order 
globule-coil transition. The degree of ionization (Fig. 3] 
(b)) follows a similar trend with the jump discontinuities 
appearing at the same value of / at which the globule- 
coil transition occurs. We note that the trends of li and 
a vs. / are the same as those observed in simulations. A 
comparison of the terms in the total free energy suggests 
that a is coupled to h through the electrostatic inter- 
action term in Eq. [5) For w < 0, the two body term 
in Eq. [5] tends to minimize h, and competes against the 
stretching term in Eq. 1111 When / > fc, h increases 
abruptly to a high value and the electrostatic interaction 
term in Eq. [5] becomes negligible. Then, a is determined 
by a competition between Fads and mobile counterion 
entropy. 

In conclusion, we have explored the globule-coil tran- 
sition in polyelectrolytes under mechanical tension. We 
have shown using Langevin dynamics simulations that 
at a critical value of the applied force, the globule-coil 
lengthening transition of the chain is accompanied by a 
correlated counterion desorption. In addition, our the- 
oretical calculations indicate that the jump in the ion- 
ization of the chain at the globule-coil transition arises 
due to a gain in mobile counterion entropy that is greater 
than the chain-counterion adsorption energy. Our results 
can help understand conformational changes in biopoly- 
mers where electrostatic effects are important, and also 
aid interpretation of single molecule experiments involv- 
ing polyelectrolytes. 
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